Veze izmedju elemenata
clear
variables
syms
omega L
C
w
w = omega
w =
ω
C = 1/(L*w^2)
C =
Pobuda
syms
Um s t
Ug = laplace(Um*sin(w*t), t, s)
Ug =
Jednacine
syms
I0 IL U0 UC
jednacine = [Ug == L*(s*IL - I0) + UC, IL == C*(-U0 + s*UC)]
jednacine =
odziv = solve(jednacine, [UC, IL])
odziv =
struct with fields:
UC: (I0*L*omega^4 + Um*omega^3 + I0*L*omega^2*s^2 + U0*omega^2*s + U0*s^3)/(omega^2 + s^2)^2 IL: (I0*L*omega^2*s - U0*omega^2 + Um*omega*s + I0*L*s^3 - U0*s^2)/(L*(omega^2 + s^2)^2)
odziv.UC
ans =
Inverzna Laplasova transformacija
syms
uC(t)
uC(t) = ilaplace(odziv.UC, s, t)
uC(t) =
vrednosti = [I0 == 0.2, L == 0.1, U0 == 2, Um == 2, w == 2*pi*50]
vrednosti =
uC_num(t) = subs(uC(t), lhs(vrednosti), rhs(vrednosti))
uC_num(t) =
Crtanje grafika
fplot(t, uC_num(t), [0, 0.7])
xlabel(
't'
);
ylabel(
'uC(t)'
);
grid
on
;